COMPUTATION OF ANALYTIC CAPACITY AND 
APPLICATIONS TO THE SUBADDITIVITY PROBLEM 



MALIK YOUNSI AND THOMAS RANSFORD 

Abstract. We develop a least-squares method for computing the analytic 
capacity of compact plane sets with piecewise-analytic boundary. The method 
furnishes rigorous upper and lower bounds which converge to the true value 
of the capacity. Several illustrative examples are presented. We are led to 
formulate a conjecture which, if true, would imply that analytic capacity is 
subadditive. The conjecture is proved in a special case. 



1. Introduction 

Let K be a compact subset of C and let f2 be the complement of K in the Riemann 
sphere, i.e. il := Coo \ K. The analytic capacity of K is 

j(K) := sup{|/'MI : / e ||/||oo < !}• 

Here /'(oo) := lim^oo z(f(z) — /(oo)) denotes the coefficient of z^ 1 in the Lau- 
rent expansion of f(z) near infinity, and H°°(£l) denotes the class of all bounded 
holomorphic functions in Q. 

Analytic capacity of compact sets was first introduced by Ahlfors pQ, in order to 
study Painleve's problem of finding a geometric characterization of the compact sets 
K that have the property that every bounded holomorphic function in Q, is constant. 
These compact sets are called removable and are precisely those of zero analytic 
capacity. Painleve's problem has been extensively studied in the last decades and is 
now considered solved. See e.g. [T7] for a survey of Painleve's problem and related 
results. 

The study of analytic capacity became even more interesting with Vitushkin's work 
on uniform rational approximation |20j . Vitushkin showed that analytic capacity 
plays a central role in the theory of uniform rational approximation of holomorphic 
functions on compact subsets of the plane. See e.g. [35] for a survey of the appli- 
cations of analytic capacity to this type of problem. 

We also mention that analytic capacity is used in fluid dynamics to study the 2- 
dimensional velocity fields induced by several obstacles, see e.g. [ID] . 
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In this article, we are primarily interested in the computation of analytic capacity. 
More precisely, our main objectives are 

• To obtain a quick and efficient method to compute the analytic capacity of 
"nice" compact sets. 

• To use this method to investigate the (still open) subadditivity problem for 
analytic capacity. 

The article is structured as follows : 

Section [2] contains the necessary preliminaries on analytic capacity. In Section[3l we 
obtain some estimates for the analytic capacity of a compact set with C°° boundary. 
Then, in Section 01 we prove that the same estimates remain valid in the case of 
compact sets with piecewise-analytic boundary (subject to suitable modifications). 
The proof relies on properties of the Smirnov classes E p (£l) on finitely connected 
domains. In Section [5l we use these estimates to obtain a numerical method for 
the computation of analytic capacity. The method gives upper and lower bounds 
for the analytic capacity j{K), and we prove that these bounds can be made arbi- 
trarily close, thus converging to the true value of j(K). In Section [6J we present 
several numerical examples, in the case of analytic boundary as well as in the case 
of piecewise-analytic boundary. The last two sections, Section [7J and Section [5J 
are dedicated to the study of the subadditivity problem for analytic capacity. The 
problem is the following : 

Is it true that 7 is subadditive, in the sense that 

1 {E\JF) < j(E)+<y(F) 
for all compact sets E,F C C ? 

Vitushkin conjectured in [20] that analytic capacity is semi- additive, i.e. 

-f(EUF) < C(j(E)+j(F)) 

for some universal constant C. He gave various applications of this inequality to 
rational approximation. Semi-additivity of analytic capacity was proved by Tolsa 
[15] , In fact, Tolsa proved that analytic capacity is countably semi-additive. How- 
ever, it is still unknown whether or not one can take the constant C equal to I . 
Davie [4] gives some applications of the subadditivity of analytic capacity to ratio- 
nal approximation theory 

In Section [71 we first use a discrete approach to analytic capacity (introduced by 
Melnikov [9]) to prove that the subadditivity of analytic capacity is equivalent to 
the subadditivity in the special case where the compact sets E, F are disjoint finite 
unions of disjoint disks, all with the same radius. This result is quite convenient 
because the numerical method described in Section [5] is very efficient for computing 
the analytic capacity of such compact sets. Then we use a discrete version of 
analytic capacity, also introduced by Melnikov 9 , to obtain a result regarding the 
behavior of the ratio 

7(£Uf) 
1 (E)+ 1 (F) 

as r — > 0, where E and F are disjoint finite unions of disjoint disks, all with radius r. 
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Finally, in Section |8l we formulate a conjecture based, among other things, on nu- 
merical evidence. A proof of this conjecture would imply that analytic capacity is 
subadditive. We end the section by giving a proof in a special case. 

In the article, we shall use the letter SI to denote a domain in the Riemann sphere, 
that is, an open and connected subset of Coo. Furthermore, when SI is said to be a 
finitely connected domain with analytic (respectively piecewise-analytic) boundary, 
we mean that the boundary of SI consists of a finite number of pairwise disjoint 
analytic (respectively piecewise-analytic) Jordan curves. Finally, we shall use A(Q) 
to denote the set of complex- valued functions that are continuous in SI, the closure 
of SI in Coo, and holomorphic in SI. 

2. Preliminaries on Analytic Capacity 

In a sense, analytic capacity measures the size of a set as a non-removable singularity 
for bounded holomorphic functions. A direct consequence of the definition is that 
7 is monotone : 

K x C K 2 => j(Ki) < j{K 2 ). 
It is also easy to prove that 

j(aK + b) = \a\j{K) 
for every a, b €E C and compact set K. In particular, 7 is invariant under translation. 

Analytic capacity is also outer regular, in the sense that, if 

K x 2 K 2 D K 3 D . . . 

is a decreasing sequence of compact sets, and if K := C\ n K n , then 7(-Kn) — > l{K) 
as n — > 00. 

It is well known that, for every compact set K, there exists an extremal func- 
tion / for j(K), that is, a function / holomorphic in SI with |/| < 1 in SI and 
/' (00) = j(K). In the case 7 (if) > 0, this function / is unique in the unbounded 
component of SI and is called the Ahlfors function for K. One verifies easily that 
the Ahlfors function / vanishes at 00. 

From Schwarz's lemma, it follows that, if K is connected, then / is the conformal 
map of SI onto the unit disk P with /(oo) = and /'(oo) > 0. As a consequence, 
we get that the analytic capacity of a closed disk equals the radius, and the analytic 
capacity of a closed line segment equals a quarter of the length. See e.g. [8]. 

2.1. Finitely connected domains with analytic boundary. Let if be a com- 
pact set in the plane and again denote by SI the complement of K in Coo- The 
following theorem says that, under certain assumptions on K, the Ahlfors function 
behaves nicely : 

Theorem 2.1 (Ahlfors 1 ]). Let us assume that SI is a finitely connected domain 
whose boundary consists of n Jordan curves. In this case, the Ahlfors function f is 
an n-to-1 branched covering of S7 onto the unit disk. Moreover, 

(1) / extends continuously to dfl, the boundary o/Sl, so that f € -A (SI); 

(2) l/l = 1 on dQ; 
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(3) / maps each of the n boundary curves homeomorphically onto the unit 
circle. 

By the Schwarz reflection principle, if in addition each boundary curve is analytic, 
then / extends analytically across the boundary. 



One way to prove Theorem 1 2. II is to use the following result : 

Theorem 2.2. Suppose that £1 is a finitely connected domain with analytic bound- 
ary. Then there exists a holomorphic function ip in CI which is the unique solution 
to the dual extremal problem 

1 



\iP(0\\d(\=inii \h(()\\dC\:heA({l),h(oc)= . 
on [J on 27ri 

Moreover, ip has the following properties : 

(1) ip G A(Q) and extends analytically across dCl; 

(2) V(oo) = 1/2™; 

(3) ip represents evaluation of the derivative at oo, in the sense that, for all 

geA(Ci), 

«/(oo) = f g((MOd(; 

JdQ. 

(4) s dn m)\m = i(K); 

(5) The extension of ip has an analytic logarithm. In particular, there exists a 
function q € A(Cl) such that q(co) = 1 and 

q(zf = 2mtf)(z) (z e 0). 

The above is essentially due to Garabedian [7 . The function rp is usually called the 
Garabedian function for f2. See also [U Theorem 4.1]. 

We end this subsection by remarking that the function q in the above theorem is, 
up to a multiplicative constant, a reproducing kernel for the Hilbert space H 2 (il). 

Indeed, recall that, for any domain ft and for < p < oo, the Hardy space H p (£l) 
is the class of all functions h holomorphic in such that the subharmonic function 
\h\ p has a harmonic majorant. This definition is conformally invariant and coincides 
with the classical one when Q is the unit disk. If Q is a finitely connected domain 
with analytic boundary and if oo £ Q, then H 2 (Q) is a Hilbert space, with respect 
to the scalar product 



(g, h ) = / g{z)h{z)\dz\, 
J on 

in which evaluation at oo is continuous. Hence there is a unique function S(z, oo) G 
H 2 (fl), called the Szego kernel function for oo, such that 



fl (oo)= / g(z)S(z,oo)\dz\ (g € H (CI)). 
Jdn 



If q is the function in (5) of the above theorem, then we have 

q(z) = 2irj{K)S(z, oo) (z G CI). 
See e.g. (HJ Theorem 4.3]. 
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2.2. Transformation of the Ahlfors and Garabedian functions under con- 
formal mapping. The goal of this subsection is to describe how the Ahlfors and 
Garabedian functions transform under conformal mapping. 

As before, let K be a, compact set in the plane, and set Q :— Coo \ K. Suppose 
that is a finitely connected domain whose boundary consists of a finite number 
of pairwise disjoint Jordan curves. 

By repeated applications of the Riemann mapping theorem, there exists another 
compact set K whose complement 17 is a finitely connected domain with analytic 
boundary and is conformally equivalent to ft. Denote by F : ft — > f2 the conformal 
map thereby obtained, normalized so that ^(oo) =00. 

It is well known that every conformal map of a Jordan domain onto the unit disk 
extends to a homeomorphism of the closure of the domain onto the closed unit disk. 
Hence, by construction, F extends to an homeomorphism of Q onto Q. Write 

1 . a_i a_ 2 
F(z) = aiz + a H 1 5- + . • . 

z z z 

near infinity. The following proposition relates the Ahlfors functions for K and K: 



Proposition 2.3. Let /, / be the Ahlfors functions for K and K respectively. Then 
j(K) = j(K)/\ai\ and the following diagram commutes, up to a multiplicative con- 
stant of modulus 1 : 



n — ^ n 



Proof. Note that / o F e H°°(n) and ||/ o F^ < 1. Thus, 

|(/oF)'(oo)| < 7 (Jf). 




However, we have 



and so 



Repeating this with / o F 1 instead gives the reverse inequality. 

By uniqueness of the Ahlfors function, we have / o F = Xf for some constant A 
with |A| = 1. 

□ 
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For the transformation of the Garabedian function ip, we need additional assump- 
tions on the boundary of f2. The reason behind this will be clear soon. 

Accordingly, we shall assume that F is C°° on the boundary. This will be the case, 
for example, if all the boundary curves of are C°° . This is a consequence of the 
following result, which dates back to Painleve's doctoral thesis : 

Theorem 2.4 (Painleve). Let D be a bounded Jordan domain with C°° boundary, 
and let f be a conformal mapping of D onto the unit disk D. Then f is C°° on D, 
the derivative f does not vanish on D, and / _1 is C°° on the closed unit disk. 

Proof. For a proof, see e.g. (31 Theorem 8.2]. □ 

We shall also need the fact that F' has an holomorphic square root in ft. This is a 
consequence of the following : 

Theorem 2.5 (Bell). Let / : f2i — > O2 be a conformal mapping between bounded 
finitely connected domains with C°° boundaries. Then f is C°° on f2i and f is 
nonvanishing on fix- Consequently, / _1 is C°° on f2 2 . Furthermore, f is equal to 
the square of a function C°° on Qx and holomorphic in Q-y. 

Proof. See [2 Theorem 12.1]. □ 
We can now prove : 

Theorem 2.6. Suppose that Q is a finitely connected domain with C°° boundary, 
and let f2 be as in the above. Let ip be the Garabedian function for f2, as in Theorem 
\2.Sl Define a function ip in fl by 

ip := —$oF)F'. 
ai 

Then ip has the following properties : 

(1) ip e A(fi); 

(2) i)(oo) = l/2wi; 

(3) ip represents evaluation of the derivative at 00, in the sense that, for all 
9 G A(fi), 

g'(oo) = [ 5 (C)^(C)dC; 

JdQ. 

(4) f dn m)\\d(\=l(K); 

(5) ip has an analytic square root in tt. More precisely, there exists a function 
q 6 A(Q,) such that q(oo) — 1 and 

q{z) 2 = 2wiip(z) (z e fi). 
Proof. (1) is clear, ip is even C°° on il. 

To prove (2), note that 

F (z) — ai 3- - . . . 

z z z- 3 

near 00, so that F'(z) — > a\ as z — > 00. Since F fixes 00, we have 

(ip o F)(oo) = ip(oo) = — 
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and (2) follows. 

To prove (3), we use the change-of- variables formula as found in e.g. [TH Theorem 
7.26]. Let g £ A(Q). Then joT 1 £ A(h) so that, by Theorem E2J 

(goF^'M = (g{F- 1 {w))i!{w)dw 
Jan 

g(z)^(F(z))F'(z)dz 

an 

= ai / g(z)ip(z)dz, 
Jan 

Here we used the change of variable w — F(z), which is legitimate since F is injec- 
tive and C°° on d£l. 



Now write 



g(z) = 5 (oo) + + ^ + 

z z z 



near 00. Then 

(g o F _1 )'(oo) = lim w(g{F~ 1 (w)) - ff(oo)) = 015(00) 

w— >oo 

and (3) follows. 

To prove (4), note that 

Mz)\\dz\ = 7^7 / |V~(^))II^(^)^I 

an l°i| Jan 

1 



_ |^(w)||gH 
l°i| Jan 



by Proposition 12.31 

Finally, (5) follows directly from Theorem 12.21 and the remark preceding Theorem 
1231 

□ 

Adopting the terminology already used in the case where fl had analytic boundary, 
from now on we shall call ip the Garabedian function for Q. 



3. Estimates for Analytic Capacity in the case of C°° boundary 

In this section, we obtain some estimates for the analytic capacity of a compact 
set whose complement is a finitely connected domain with C°° boundary. First, we 
need a lemma : 
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Lemma 3.1. Let K be a compact set in the plane whose complement ft is a finitely 
connected domain with C°° boundary. If f G -A(ft), then 



f'(oo) = 




Proof. Write 

f{z) = /(oo) + + — + ... 

z z z 

near oo. Let C be a circle centered at the origin and containing K, with radius 
sufficiently large so that the above expression for / holds on C. Then we have 



Indeed, the above is clear if / is holomorphic in a neighborhood of ft, and such 
functions are uniformly dense in A(Q,), by Mergelyan's theorem. 

On the other hand, the right-hand side in the last expression is just f'(oo). This 
can be seen by substituting the power series expression for / into the integral and 
integrating term by term. 

□ 

Now we can prove : 

Theorem 3.2. Let K be a compact set in the plane, and suppose that the comple- 
ment Q of K is a finitely connected domain with C°° boundary. Then 



(1) 7 (K) = min ji- J \g(z)\ 2 \dz\ : g G A(Sl) , g(oc) = 1 j 

and 



(2) j(K) = maxi2Re/i'(oo) - — / \h(z)\ 2 \dz\ : h € A(Q) , h(oo) = 1 . 

I 27r Jan J 

Here the minimum and maximum are attained respectively by the functions g = q 

and h = fq, where f is the Ahlfors function for K and q is the function of Theorem 

m 

Estimate ([IJ was already known in the case of analytic boundary, since the work 
of Garabedian [7]. It is usually referred to as Garabedian's duality. 

Proof. Let / be the Ahlfors function for fi, so that / e A(fl) with /(oo) = and 
|/| e 1 on dil and /'(oo) = ^(K). Let ip be the Garabedian function for il in the 
sense of Theorem 12 .61 and denote by q the function in A(Q) with q(oo) = 1 and 

q(z) 2 = 2niip(z) (z G ft). 



To prove (P, let g G A(ft) with 5(00) = 1. We have 

7(Jf) = /'(oo) = (/ 5 2 )'M = / f{z)g(z) 2 dz 

2m Jon 
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by Lemma [3TT1 Thus, 

l(K) < I- [ |/(z)||g(z)| 2 |^l = ^ / \9{z)\ 2 \dz\ 
2tt Jon Ztt Jgn 

Taking the minimum over all such functions g, we obtain 



±- Jjg(z)\ 2 \dz\ :g£ A(Q), ff (oo) = l| 



(3) 7 (if)<min |) 

On the other hand, take g :— q. Then g £ A(Q) with g(oo) — 1 and 

±[ \ g (z)\ 2 \dz\ = ± [ \q{z)\ 2 \dz\ = / |#*)IH=7(*"), 
Z7r Jan Z7r Jan Jan 

by Theorem 12.61 Combining this with inequality we obtain (fl}. 



For ([2]), consider the function h = fq. Then h £ A(fl), h(oo) = and 

27 



2Re//(co)-^- / \h{z)\ 2 \dz\ = 2 1 {K)- 1 (K)= 1 {K). 
'an 



Thus 



max < 2 Re h'(oo) / 

I V ' ^ hi 



(4) max<{ 2Re/i'(oo) - — / |/i(z)| : h £ A(Q) , h(oo) = > > j(K). 

'' I an J 

On the other hand, let h £ A(il), h(oo) = 0. Denote by T(z) the unit tangent 

vector to dVL at z, that is, dz = T(z)\dz\ with |T| = 1. 

Let (hi,h 2 ) denote J dn hi(z)h 2 (z)\dz\ and ||/i||2 := (h 7 h). 

Then 

<^- \\h- iqT\\l= ±-\\m + ^\\q\\i + 2Re±-{h,-iqT), 

Z7T Z7T Z7T Z7T 

so that 

0<-^/ |Mz)| 2 |dz|+ 7 (^)-2Re-^ / h(z)q(z)dz. 
27r Jan 2tti ,/ an 

By Lemma [5TT1 it follows that 

0<— / \h(z)\ 2 \dz\+j(K)~2Re{h q y{oo). 
2tt Jao 

Since (?(oo) = 1 and /i(oo) = 0, we have (/ig)'(oo) = ft'(oo), and thus 

j{K)>2Reh'(oo)- — I \h(z)\ 2 \dz\. 

2tt Jan 

Combining this with inequality we obtain (|2|). 

□ 

4. ESTIMATES FOR ANALYTIC CAPACITY IN THE CASE OF PIECEWISE- ANALYTIC 

BOUNDARY 



The objective of this section is to extend the estimates of Theorem 13.21 to another 
interesting case, that of sets with piecewise-analytic boundaries. 

Let us assume that K is a compact set such that f2 :— Coo \ if is a finitely connected 
domain with piecewise-analytic boundary. By this, we mean that the boundary con- 
sists of a finite number of non-intersecting Jordan curves, and that each boundary 
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curve is the union of a finite number of analytic arcs. We further assume that every 
intersecting pair of analytic arcs meets at a corner that is conformally equivalent to 
a sector. More precisely, we suppose that if two analytic arcs intersect at a point w, 
then there exists a conformal map defined in a neighborhood V of w, and mapping 
V fl fi onto a sector {re ie : < r < 1, < 9 < a}, where < a < 2ir. 

To obtain the estimates of Theorem 13. 2\ we used the fact that the functions q and 
fq extend continuously to the boundary, where / is the Ahlfors function for K and 
q is the square root of 2m times the Garabedian function ip. 

In the case of piecewise-analytic boundary, this remains true for the Ahlfors function 
/. However, there are some issues regarding the Garabedian function : tp will have 
discontinuities at the finite set E where the boundary fails to be smooth [E is made 
of the endpoints of the analytic arcs in the boundary) . Consequently, in order to 
extend Theorem 13 . 21 we need to replace A(fi) by a larger class of holomorphic 
functions in f2. It turns out that Smirnov classes over finitely connected domains 
are precisely what we need. 

4.1. Smirnov classes on finitely connected domains. Our objective now is to 
present the theory of Smirnov classes E p (£l) on finitely connected domains. For 
more details, we refer the reader to j6] and [T3] , 

Let Q be a finitely connected domain with rectifiable boundary. By that, we mean 
that dfl consists of a finite number of pairwise disjoint rectifiable Jordan curves. 
Let 1 < p < oo. 

We say that a function h belongs to the Smirnov class E p (£l) if h is holomorphic in 
f2 and if there exists a sequence {£l n } of finitely connected subdomains of with 
rectifiable boundaries {C n } such that : 

(1) f2„ eventually contains each compact subset of Q, 

(2) limsup^^ J Cn \h{z)\ p \dz\ < oo, 

(3) the lengths of the curves of the C n 's are uniformly bounded. 

In the simply connected case, it is well known that condition (3) is a superfluous 
requirement in the definition. This remains true in the finitely connected case. 

It is not hard to prove that A(f2) C E p (Cl). Moreover, it is also well known that 
E p (£l) reduces to the classical Hardy space H p (il) when fl is the unit disk, and this 
is also true if SI is a finitely connected domain with analytic boundary. However, in 
general, neither of the inclusions hold. Even for simple cases like Jordan domains 
with polygonal boundaries, the two classes are not equal. 

The following well-known result is a generalization of Fatou's theorem on the clas- 
sical Hardy spaces. 

Theorem 4.1. Let p > 1 and let fl be a bounded finitely connected domain with 
rectifiable boundary. Suppose that h € E p (fl). Then 

(1) h has nontangential boundary values h* almost everywhere on dQ, and 
h* G L p {dn). 
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(2) h is the Cauchy integral of h* : 

Hz) = ±- t p^dc (z e n). 

2tt« Jan C - z 

We shall also need the following generalization of Lemma 13.11 : 

Corollary 4.2. Let K be a compact set in the plane, and suppose that the com- 
plement of K is a finitely connected domain with rectifiable boundary. Let 
h € E^fl). Then 

h'(oo) = -L / h*(()d(. 
2m Jan 

Proof. Translating K , we may suppose that it contains in its interior. Set fio := 
{z : l/z G f2}. Then SIq is a bounded finitely connected domain with rectifiable 
boundary and G Qq. Define ho in f^o by ho(z) — h(l/z). It is easy to see that 
ho € E 1 (flo), since if To is any rectifiable Jordan curve in Qo not passing through 

0, then 

/ \h (z)\\dz\ = / IM«0lS' 

Jr Jr \ w \ 

where T := {z : l/z G r }, and the function 1/w is bounded in ft. 

Now, the function go(z) := (ho(z) — ho(0)) / z clearly also belongs to i? 1 (i7o), so that 

h (oo) = go{0) = — / dz = — / — —dz, 

2tt« JdQ z 2m Jdn z 

where we used the preceding theorem. Making the change of variable w — l/z, we 
obtain the result. 

□ 

4.2. The Garabedian function. Let us now return to the case where K is a 
compact set in the plane whose complement 17 is a finitely connected domain 
with piecewise-analytic boundary. In this case, there are some issues regarding 
the Garabedian function ip. If we proceed as in Theorem 12.61 and define 

ib := —(^oF)F', 
ai 

then i/j will not extend continuously to the boundary : the first factor ip o F is in 
A(iT), but the second one F 1 has singularities at the endpoints of the analytic arcs, 

1. e. at the points of E. However, F extends analytically across any analytic arc in 
the boundary, so ip as defined is continuous in Q \ E. 

The following result is the analogue of Theorem l2.6l in this new setting. Recall that 
n is a finitely connected domain with analytic boundary conformally equivalent to 
n, and that F : Q — > is a conformal map, normalized so that F (oo) = oo, and 
with expansion 

b [z] = a x z + a -\ 1 5- + . . . 

z z z 

near infinity. 

Theorem 4.3. Let if be the Garabedian function for Q. Define a function ip in Q 
by 

xj; : = —U)oF)F'. 
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Then ip is a Garabedian function for SI, in the sense that 

(1) ip is holomorphic in SI and continuous in Sl\ E, 

(2) i){oo) = l/2m, 

(3) ip represents evaluation of the derivative at oo, in the sense that for all 
9 € A(Q), 

ff'(oo) = f ,g(CMCR, 
Jon 

(4) f aa MC)\\dC\=?(K), 

(5) ip has an analytic square root in SI. More precisely, there exists a function 
q holomorphic in SI and continuous in Sl \ E such that q(oo) — 1 and 

q{zf = 2tt#(X) (zeU\E). 

(6) ip € E 1 (Sl). Consequently, q € S 2 (Q). 

For the proof, we need an analogue of Theorem 12.51 : 

Lemma 4.4. Let Q, f2 and F : Q CI be as in the above. Then F' has an analytic 
square root in SI. More precisely, there exists a function h holomorphic in SI and 
continuous in Sl \ E such that 

h(zf = F'(z) (z€Q\E). 

Proof. The proof of [3l Theorem 12.1] also works in our case. However, we can use 
the fact that F is a composition of Riemann maps to obtain a more elementary 
proof, as follows. 

Let n be the number of curves in the boundary of SI. Recall that by construction, 
F is a composition of n Riemann maps : 

F = 4>\ o 02 o • • • o 4> n , 

where each (pj maps some unbounded Jordan domain onto C 00 \ID), with </>j(oo) = oo. 
We proceed by induction on n. 

First, consider the case n = 1. Translating SI if necessary, we can suppose that 
£ Si. Put D x :={z:l/ze SI} and define 

G{Z) --=WJZ) iZGDl) - 
Then G is a conformal mapping of the bounded Jordan domain Di onto ID with 
G(0) = 0. We know that G extends to a homeomorphism of D\ onto D and 
analytically across any analytic arc of the boundary. Note that the boundary of D\ 
consists of a finite number of analytic arcs separated by a finite set of points. Call 
this finite set of points E. It is easy to see that G extends analytically to a simply 
connected domain U containing D\ \ E, with G' ^ there. 
Write G'(z) = g(z) 2 for some g holomorphic in U. We have 



so that 



: '"■■' 3 
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But both F' and h are holomorphic in a neighborhood of fl \ E. In particular, we 
have 

F'{z) = h{zf (zeTi\E). 

This completes the proof for the case n = 1. 

Now, suppose that the result holds for n—1, where n > 2. Since F — 0io0 2 °- • -°0 n , 
we have 

F'(z) = 4>[((4> 2 o • • • o 4> n ){z)){<}> 2 o . . . 4> n )'(z) 

and the result follows from the induction hypothesis and the case n = 1. 

□ 

We can now proceed to the proof of Theorem 14.31 : 

Proof. We already know that (1) holds, and the proof of (2) is exactly the same 
than the one in Theorem 12.61 Moreover, since F is differentiable everywhere on 
the boundary except at a finite set of points, we can use the change-of-variablcs 
formula found e.g. in [141 Theorem 7.26]. Points (3) and (4) then follow exactly as 
in the proof of Theorem 12.61 

Point (5) follows directly from Theorem 12. 21 together with Lemma T4.4I 

For (6), since if> := F)F' and ip o F is bounded in f2, it suffices to show that 
F' € But this is clear, since if C is any rectifiable Jordan curve in fi, then 

/ \F'(z)\\dz\ = [ \dw\. 

JC JF{C) 

□ 

4.3. Proof of the estimates. We can now prove : 

Theorem 4.5. Let K be a compact set in the plane, and suppose that the com- 
plement £1 of K is a finitely connected domain with piecewise- analytic boundary. 
Then 



j(K) = min 




and 



j(K) = max {2 Re h'(oo) / \h*(z)\ 2 \dz\ : h G E 2 {Vl) , h(oo) = 

I 27r Jan 

Here the minimum and maximum are attained respectively by the functions g = q 
and h = fq, where f is the Ahlfors function for K and q is the function of Theorem 

m 

Proof. The proof is exactly the same as the one in Theorem 13.21 Use Theorem 14.31 
instead of Theorem 12.61 and Corollary 14.21 instead of Lemma 13.11 

□ 
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Remark 1. It follows from the first estimate that the function q is unique and, 
consequently, the Garabcdian function tp too. Indeed, q is an clement of minimal 
norm in the convex set S := {g G E 2 (fl) : g(oo) — 1}, which is necessarily unique 
by an elementary Hilbert-space argument. 

5. Computation of Analytic Capacity 

5.1. Description of the method. In this section, we present a method based on 
the estimates of Theorem 13.21 (respectively Theorem I4.5[) to compute the analytic 
capacity of a compact set K whose complement is a finitely connected domain 
with analytic (respectively piecewise-analytic) boundary. The method yields upper 
and lower bounds for j{K). 

Let Ao(O) := {/ e A(Cl) : f(oo) — 0} and let JF be a subset of Aq(Q) whose span is 
dense, with respect to the L 2 -norm on d£l. For example, T could be the set of all 
functions of the form (z — a)~ n , where n <E N and a belongs to some prescribed set S 
containing one point in each component of the interior of K . This is a consequence 
of Mergelyan's theorem. 

Let A = {gi, <?2; • • • > <7n} be a finite subset of J- . The functions g\, g%, . . . , g n will be 
called approximating functions. 

The method for the upper bound is based on the following : 

• Finding the function g in the span of gi, g%, . . . , g n that minimizes the quan- 
tity 

i- / |l+.g(z)| 2 |dz|. 
i7r Jon 

In view of Theorems 13.21 and 14.51 this gives an upper bound for 7 (if). 
More precisely, the method is the following. 

• Define 

n 

i=i 

where a\, . . . , a n are complex numbers to determine. Write ctj := Cj + idj, 
and then compute the integral 

i / |l+ 5 (z)| 2 |dz|. 
2lT Jon 

This gives an expression that can be written in the form |x T Ax + bx + c, 
where x = (ci, c%, ■ . . , Cn, d%, d,2, ■ ■ • , d n ), and where A is a real symmetric 
positive-definite 2n x In matrix, b is a real vector of length 2n and c is a 
positive constant. 

• Find the Cj 's and dj 's that minimize this expression. This can be done for 
example by solving the linear system 

Ax + b = 0. 

• Create a new set of approximating functions A by adding functions from 
T to A, and then repeat the procedure with A replaced by A. 
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The above yields a sequence of decreasing upper bounds for 7 (if). 

Clearly, it can be adapted to yield a sequence of increasing lower bounds for j{K), 
using the other estimate of Theorems 13.21 and 14.51 

5.2. Convergence of the method. In this subsection, we prove that the upper 
and lower bounds obtained with the method can in principle be made arbitrarily 
close. 

First, recall that the minimum and maximum in Theorems 13.21 and 14.51 are attained 
respectively by the functions q and fq, where / is the Ahlfors function for K and 
q is the square root of 2iri times the Garabedian function for fl. 

If the boundary of £1 is C°°, then both of these functions belong to A(SY). It 
follows from Mergelyan's theorem that we can approximate them uniformly on d£l 
by rational functions with poles in S 1 , where S is some prescribed set containing at 
least one point in each component of the interior of K. On the other hand, Lemma 
O implies that if h,R € A(Q) and \h- R\ < e on dtt, then \R'(oo) - h'(oo)\ < Ce 
where C depends only on Q. 

Hence, it follows that for every e > 0, there exist rational functions R\, R2 vanishing 
at 00 and with poles in the prescribed set 5, such that 



This proves the convergence of the bounds in the C°° boundary case. 

For the piecewise-analytic boundary case, we need an analogue of Mergelyan's The- 
orem for the Smirnov class E 2 (Vl). 

Assume that if is a compact set in the plane whose complement is a finitely con- 
nected domain with piecewise-analytic boundary. 

We know that E 2 (Q) contains A(fl), but is it true that A(il) is dense in E 2 (p)l 
In other words, can every function h in E 2 (fl) be approximated on the boundary 
by functions /„ in A(Cl), in the sense that 



as n — > 00? 

In turns out that the answer is yes. Before we prove this, we need the definition of 
Smirnov domains : 

Let D C C be a bounded Jordan domain with rectifiable boundary. Since D is 
simply connected, there is a conformal mapping <f> of the open unit disk D onto D. 
It is well known that <f>' is in ii^O), and, since it has no zeros, we have a canonical 
factorization of the form 




1 



[ \1 + R 2 {z)\ 2 \dz\~ e < -f{K) < 2Rc R[(w) - ^- [ \R X {z)\ 2 \dz\ 



+ e. 




<P'(z) = S(z)Q(z) (z e ©), 
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where S is a singular inner function and Q is outer. We say that D is a Smirnov 
domain if S = 1, that is, if <j)' is outer. It can be shown that this definition is 
independent of the function cf>; it depends only on the domain D. A simple suffi- 
cient condition for D to be a Smirnov domain is that arg <j)' be bounded either from 
above or below. Geometrically, this means that the local rotation of the mapping 
is bounded; loosely speaking, the boundary curve cannot spiral too much. In par- 
ticular, D is a Smirnov domain if it has smooth (or piecewise-smooth) boundary. 
We refer the reader to [5] or [TTJ Chapter 7] for more details on Smirnov domains. 

We say that a function h £ L p (dD) belongs to the L p (9D)-closure of the polyno- 
mials if there is a sequence (p n ) of polynomials such that 



It is convenient to identify E P (D) with its set of boundary values functions. Thus, 
E P (D) is a closed subspace of L p (dD) which contains the polynomials, hence also 
their closure. For the reverse inclusion, we have the following criterion: 

Theorem 5.1. Let D be a bounded Jordan domain with rectifiable boundary, and 
let 1 < p < oo. Then E P (D) coincides with the L p {dD)- closure of the polynomials 
if and only if D is a Smirnov domain. 

Proof. See e.g. [6J Theorem 10.6]. □ 

Our objective is to generalize Theorem 15. II to finitely connected domains. We shall 
need the following theorem : 

Theorem 5.2 (Decomposition Theorem). Suppose that D is a bounded finitely 
connected domain whose boundary consists of pairwise disjoint rectifiable Jordan 
curves T\, I^, . . . , r n , where the outer boundary of D is T±. For 1 < j < n, let Dj 
be the component of Coo \ Fj that contains D. 

Let 1 < p < oo, and let h £ E P (D). Then h can be decomposed uniquely as 



The following is a generalization of Theorem 15. II to finitely connected domains, in 
the case p = 2 : 

Theorem 5.3. Let D and Dj (1 < j < n) be as in Theorem \5.2\ and suppose in 
addition that the curves Tj are piecewise analytic. Let a\ '■= oo and for 2 < j < n, 
fix a point aj in the interior of the complement of Dj . 

Then the rational functions with poles in the set {a±, 02, ... , a„} are dense in E 2 (D). 
In other words, for every h £ E 2 (D), there exists a sequence (R n ) of rational 
functions with poles in the prescribed set {a±, 02, . . . , a n } such that 




(n — > 00). 



h{z) = h\(z) + h 2 (z) H \-h n (z) (z £ D) 

where each hj belongs to E p (Dj) and hj (00) = for 2 < j < n. 

Proof. See [19] . 



□ 




as n — > 00. 



For the proof, we need the following lemma 
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Lemma 5.4. Let 1 < p < oo, let U be a bounded Jordan domain with rectifi- 
able boundary and let K be a compact subset of U . Then there is a constant M , 
depending only on p and K , such that 

\9{w)\ < M\\g*\\ p 
for all w G K and for every function g G E P (U), where 

n-= (J gu \g*(z)\ p \dz\y p '. 

Proof. This is a simple application of Theorem 14. II and Holder's inequality. 

□ 

We can now prove the theorem : 

Proof of Theorem ] 5. 3[ By the decomposition theorem, it suffices to show that, if 
h € E 2 (Dj), then there exists a rational function Rj with poles only at <Zj, such 
that the integral 

/ \h*(z)-R j (z)\ 2 \dz\ 

J 3D 

can be made arbitrarily small. Let e > 0. 

First consider the case j — 1, so that h G E 2 (Di). Since D\ is a bounded Smirnov 
domain with boundary Ti, we can apply Theorem 15. f I and find a polynomial P\ 
such that 

/ \h*(z)-P 1 (z)\ 2 \dz\<e. 

Now, since the function (h — Pi) G E 2 (Di), we know by Lemma [5.41 that there 
exists a constant M, depending only on the curves T3, . . . , r„, such that 



\h(w) - Pi(w)\ 2 < M / \h*{z) - Pi{z)\ z \dz\ <Me (w e T 2 U • • ■ U T n ). 
Thus, we have 

/ \h*(z)-P 1 (z)\ 2 \dz\ = I \h*(z)-Pi(z)\ 2 \dz\+f2 I \h{z)-P 1 {z)\ 2 \dz\ 

JdD JVi j=2 JTj 

< e(l+ML), 



where L is the sum of the lengths of the curves r 2 , r 3 , . . . , T n . Since the right side 
can be made arbitrarily small, this concludes the proof for the case j = 1. 



Suppose now that 2 < j < n and let h G E 2 (Dj). We can suppose h(oo) = 0, since 
it is part of the conclusion in the decomposition theorem. Translating Dj , we may 
suppose that a 3 = 0. Let Dj := {z : 1/z G Dj}. Then Dj is a bounded Jordan 
domain with piecewise-analytic boundary; in particular it is a Smirnov domain. 
Define a function g in Dj by g(z) := h(l/z). It is easy to check that g G E 2 (Dj). 
Also, g vanishes at 0, so the function g{z)/z also belongs to E 2 (Dj). We have, by 
Theorem 15.11 

9*{z) 



3-D, 



P(Z) 



\dz\ < e 
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for some polynomial P. Making the change of variable w — 1/z, we obtain 

2 



dDi 



wh*(w)-P [ — 
w 



T^\ dw \ <€ - 



Let Q(w) :— (l/w)P(l/w), so that Q is a rational function with poles only at 0. 
The above inequality becomes 

\h*(w)-Q(w)\ 2 \dw\ <€. 

Now, apply Lemma HT^l once again to conclude that there is a constant M depending 
only on the curves for k ^ j, such that 



g[z) 



P(z) 



< Me 



for all z G T/j, k ^ j, where ]?/. := {w : 1/w € Tj}. Thus, 

/ \h*(w)~Q(w)\ 2 \dw\ = f \h*{w)-Q(w)\ 2 \dw\+ V f \h(w) - Q(w)\ 2 \dw\ 

JdD JTj kltj Jr k 



< e 



E L 



P(z) 



dz 



< e + ML'e = (1 + ML')e, 



where L' is the sum of the lengths of the curves Tfc, k ^ j . Since the right side can 
be made arbitrarily small, this completes the proof of the theorem. 

□ 

Remark 2. In the above proof, piecewise-analyticity of the boundary is only as- 
sumed so that the domains Di and Dj (2 < j < n) are Smirnov domains. The 
result therefore remains true under this weaker assumption. 



Now, let us apply Theorem 15.31 to our case : i.e. K is a compact set in the 
plane whose complement fl is a finitely connected domain with piecewise-analytic 
boundary consisting of n curves. Let D be the bounded domain obtained by adding 
an outer boundary circle, say L , with radius sufficiently large so that K is contained 
in the interior D of r . If h £ E 2 (Q), then clearly h € E 2 (D). The decomposition 
theorem gives 

h = h + h\ H + h n 

where each hj belongs to E 2 (Dj). Thus, 

h Q = h-h 1 h n , 

and, since h is holomorphic in the complement of K, this gives an analytic extension 
of ho to the entire plane, bounded near oo. By Liouville's theorem, ho is constant. 
Now, since D is a bounded finitely connected domain with piecewise-analytic bound- 
ary, we can apply Theorem 15.31 to h on D. Since ho is constant, we can omit the 
pole at oo. Thus, we have proved : 
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Corollary 5.5. Suppose that K is a compact set in the plane whose complement 
51 is a finitely connected domain with piecewise- analytic boundary consisting of n 
curves. For each 1 < j < n. fix a point aj in the interior of each component of 
K . Then the rational functions with poles in the prescribed set {ai, a 2 , ■ ■ • , a n } are 
dense in E (Q). 

6. Numerical examples 

In this section, we present several numerical examples to illustrate the method. All 
the numerical work was done with MATLAB. 

6.1. Analytic boundary. 

Example 6.1.1. Union of two disks. 

Here K is the union of two disks of radius 1 centered at —2 and 2. 





0.5 - 




8 -2 -i 


1 

-0.5 - 

_ i - 





Figure 1 . The compact set K for Example 16.1.11 

A natural choice here for the approximating functions g^s is to take powers of 
l/(z — 2) and l/(z + 2). However, we shall instead consider functions of the form 



where the cij's are distinct points in the interior of K. The reason behind this is 
purely numerical : with these functions, the integrals involved in the method can be 
calculated analytically, using the residue theorem for example. This way, we avoid 
the use of numerical quadrature methods, and this results in a significant gain in 
efficiency. 

The locations of the poles aj are arbitrary. Typically, for each disk centered at c 
with radius r, we put poles at the points 

{c, c+ri, c+rii, c-ri,c-rii, c+r 2 ,c+r 2 i, c-r 2 , c-r 2 i, . . . , c+r„, c+r n i, c-r n , c-r n i], 
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where n, . . . ,r n are equally distributed between and r. 

The table below contains the bounds for j(K) obtained with the method. 

Table 1. Lower and upper bounds for j(K) for Example 16. 1.11 



Number of poles per disk 


Lower bound for 7 (if) 


Upper bound for j(K) 


Time (s) 


1 


1.875000000000000 


1.882812500000000 


0.003279 


5 


1.875593064023693 


1.875619764386366 


0.007051 


9 


1.875595017927203 


1.875595038756883 


0.012397 


13 


1.875595019096871 


1.875595019097141 


0.017422 


17 


1.875595019097112 


1.875595019097164 


0.027115 



We end this example by remarking that, in this particular case, there is a formula 
for j(K). Indeed, if 

K = B(-c,r)UB(c,r), 
where < r < c, then we have the formula 



(5) #) = ^(^"^ 2(?)2 

Here $2 is one of the so-called Jacobi theta-functions : 



Ml) :=E« (n+1/2)3=2 3 1/4 II( 1 -? 2n )( :L + 3 2n ) 2 ' 

nGZ n=l 
00 

Mq) :=xy a = Il( 1 -3 2 ")( 1 +9 2n " 1 ) 2 ' 



n£Z n=l 

00 



n&L n=l 

The argument q is given by the solution in (0, 1) of the equation 



An easy calculation gives 



2c 2 ~r 2 ~ 2CV? 
1= 



Formula ([5]) is easily deduced from a formula of Murai in [10] , by making the well- 
known change of variables 

, _ Mq) 2 

and using the identities relating theta-functions and elliptic integrals. 

(We mention though that, in the formula for 7 (if) in [10], there is a factor c missing, 

and the formula should read 
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where F is the complete elliptic integral of the first kind.) 

A simple calculation shows that formula ([5]) can also be written in the form 

j(K) = y/<?-rH 2 {q) 2 . 
Substituting c = 2 and r = 1 gives 

-y(K) « 1.8755950190971197289. 
Compare this with the bounds obtained in Table 1. 
Example 6.1.2. Union of 25 disks. 

Each disk below has a radius of 0.4. 



-2 ^ 



o o o 

o o o 
ooo ooo 

o o oo 

Q(pQ t Q t QQ 



Figure 2. The compact set K for Example 16.1.21 
Table 2. Lower and upper bounds for ^(K) for Examplc l6.1.2l 



Number of poles per disk 


Lower bound for ^(K) 


Upper bound for "f(K) 


Time (s) 


1 


4.073652478223290 


4.219704181009330 


0.177746 


5 


4.148169157685863 


4.148514554979665 


3.702191 


9 


4.148331342401185 


4.148332498165111 


11.606526 


13 


4.148331931858607 


4.148331938572625 


24.848263 


17 


4.148331934292544 


4.148331934334756 


41.342390 
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Example 6.1.3. Union of four ellipses. 

Here is another example for the computation of the analytic capacity of a compact 
set with analytic boundary. 

The compact set K is composed of four ellipses centered at —3, 3, lOi, — lOz. Each 
ellipse has a semi-major axis of 2 and a semi-minor axis of 1 : 




Figure 3. The compact set K for Example 16.1.31 
Table 3. Lower and upper bounds for j(K) for Examplc l6.1.3l 



Number of poles per ellipse 


Lower bound for 'y(K) 


Upper bound for "f(K) 


Time (s) 


1 


4.290494449193028 


5.652385361295098 


0.962078 


5 


5.252560204660928 


5.409346641724527 


17.268477 


9 


5.356419530523225 


5.377445892435984 


54.260216 


13 


5.370292494009306 


5.372648058950175 


111.424592 


17 


5.371877137036634 


5.372044462730262 


190.042871 


41 


5.371995432221965 


5.371995878776166 


1100.468881 



In this case, the integrals involved have to be calculated numerically. We used a 
recursive adaptive Simpson quadrature with an absolute error tolerance of 10~ 9 . 
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6.2. Piecewise-analytic boundary. In this subsection, we shall consider exam- 
ples of compact sets K whose complements VL are finitely connected domains with 
piecewise-analytic boundary. 

Example 6.2.1. Square. 

In this example, we consider the square with corners 1, i, — 1, —i. 




Figure 4. The compact set K for Example 16.2.11 

We fix an integer n, and then consider the approximating functions 

1 1 1 

z ' z 2 ' z n 

Below is a table containing the bounds obtained for different values of n : 



Table 4. Lower and upper bounds for *y(K) for Example 16.2.11 



n 


Lower bound for 7 (if) 


Upper bound for j(K) 


Time (s) 


2 





707106781186547 





900316316157106 


0.021981 


3 





707106781186547 





900316316157106 


0.069278 


4 





707106781186547 





887142803070031 


0.109346 


5 





746499705182962 





887142803070031 


0.145614 


6 





746499705182962 





887142803070031 


0.202309 


7 





746499705182962 





887142803070031 


0.295450 


8 





746499705182962 





881014562149127 


0.347996 


9 





761941423753061 





881014562149127 


0.414684 


10 





761941423753061 





881014562149127 


0.595552 


15 





770723484232218 





877175902241141 


2.425285 


20 





776589045256849 





872341829081944 


5.537981 


25 





784189460107018 





870656623669828 


10.002786 


30 





786857803378602 





869257904380382 


16.344379 


35 





789068961951613 





868068649269412 


26.109797 


40 





790942498354322 





866133165258689 


33.595790 
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We immediately see that the convergence is very slow, as opposed to the results 
obtained in the case of compact sets with analytic boundaries. 

The main issue here is that we do not consider the geometric nature of the boundary. 
In order to accelerate convergence, our choice of approximating functions should 
take into account the different points where the boundary fails to be smooth. 

We propose the following method : 

In view of Theorems 14.31 and 14.51 the functions that we want to approximate are 



q(z) = cyj V o F){z)^W{z) 



and 



f(z)q(z) = cf(z) X J^oF)(z)y/F 7 {z) 
for some constant c, where F is a conformal map of f2 onto Coo\D, with F(oo) = oo. 



By Theorems 12.11 and 12.21 we know that the functions / and y (t/i o F) are contin- 
uous up to the boundary. These functions can thus be approximated by rational 
functions with poles inside the square. 

All that remains is to add approximating functions that behave like y/W at the 
corners. If a is one of the corner in the boundary, then F should, in some sense, 
straighten out the angle from 37r/2 to tt, that is F(z) must behave like 

(z-a) 2 / 3 

near a. Differentiating and then taking square root, we find that (z — a) -1 / 6 should 
be, up to a multiplicative constant, a good approximation to W F'(z) near a. Since 
we want functions that are holomorphic near oo, we shall instead consider 

-1/6 



Z 

In view of all of the above, we propose the following method for the computation 
of y(K) : 

Fix an integer n. Then add the approximating functions 

IM 

z k 

for j = 0, 1, 2, 3, 4 and k = 1, 2, . . . , n, where a\, a,2, a.3, 04 are the corners of the 
square, 

fo(z) := 1 

and 

f j{z ):=(^l) V6 ( J = 1,2,3,4). 
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We now present some numerical examples to illustrate this new method. 

For the first example, we use the new method to compute the analytic capacity of 
the square of Example 16.2.11 : 



Table 5. Lower and upper bounds for j(K) 



n 


Lower bound for "f(K) 


Upper bound for j(K) 


Time (s) 


2 


0.834566926465074 


0.835066810881929 


1.334885 


3 


0.834609482283050 


0.834678782816948 


2.918624 


4 


0.834622127643984 


0.834628966618492 


5.220941 


5 


0.834626255962448 


0.834627566559480 


8.022274 


6 


0.834626584020641 


0.834627152182154 


11.542859 



Since K is connected, we have that 

j(K) = cap(if) = V2 v '' « 0.83462684167407318630, 
where cap(if) is the logarithmic capacity of K. 

Example 6.2.2. Union of two squares, one equilateral triangle and one rectangle. 




Figure 5. The compact set K for Example 16.2.21 

Our method can easily be adapted to other compact sets with piecewise-analytic 
boundary. Indeed, suppose that K is a compact set whose boundary consists of 
m piecewise-analytic curves, say 71, . . . ,7m- First, fix a point c in the interior of 
71, and let 01, 02, . . . , ajv be the different points in 71 where the curve fails to be 
smooth. Suppose that 71 makes an exterior angle of oij at the point aj, where 
< ctj < 2ir. Then we add the following approximating functions : 
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m 

(z — c) k 

for j = 0, 1, 2, . . . , N and k — 1, 2, . . . , n, where 

fo(z) ■= 1 

and 

/^) : =(7Z7j (i = l,2,...,iV). 

All that remains is to repeat the procedure for the other curves. 



Table 6. Lower and upper bounds for ^{K) for Example 16.2.21 



n 


Lower bound for j(K) 


Upper bound for "f(K) 


Time (s) 


1 


2.688593215018632 


2.724269900679792 


10.371944 


2 


2.693483826380926 


2.695819902453329 


32.881242 


3 


2.693867645864377 


2.694261483861710 


71.562216 


4 


2.693961062687599 


2.694025016036611 


122.607285 


5 


2.693971930182724 


2.693982653270314 


184.203252 



Example 6.2.3. 

Our last example is a non-polygonal compact set with piecewise-analytic bound- 
ary. It is a typical example of the kind of geometry that often arises in applied 
mathematics, featuring smoothness of the boundary with the exception of a few 
singularities. 




Figure 6. The compact set K for Example 16.2.31 

The compact set K is composed of the unit disk and two half-unit-disks centered 
at 3 and 3i. 



COMPUTATION OF ANALYTIC CAPACITY AND APPLICATIONS TO THE SUBADDITIVITY PROBLEM7 



Table 7. Lower and upper bounds for j(K) for Example 16.2.31 



n 


Lower bound for 7(F) 


Upper bound for j(K) 


Time (s) 


2 


2.118603690751346 


2.123888275897654 


2.546965 


3 


2.120521869940459 


2.121230615594293 


4.926440 


4 


2.120666182274863 


2.120803766391281 


9.488024 


5 


2.120694837101383 


2.120716977856280 


13.679742 


6 


2.120703235395670 


2.120709388805280 


22.344576 


7 


2.120704581010457 


2.120707633546616 


28.953791 


8 


2.120705081159854 


2.120706704970516 


34.781046 



The above compact set was considered in [13] and then in |12j , for the computation 
of logarithmic capacity. It was shown that cap(F) £ [2.19699,2.19881]. Our results 
are thus consistent with the well-known inequality 

j(K) < cap(-fT). 

Before leaving this section, a brief remark is in order. Comparing the results of 
Subsection 6.2 with Subsection 6.1, we see that the convergence is quite a bit 
slower in the case of piecewise-analytic boundary, compared to the case of analytic 
boundary. In fact, it is known that one cannot hope for similar convergence in 
both cases. This is related to the fact that if the boundary curves are piecewise- 
analytic but not analytic, then the extremal functions in Theorem l4.5l do not extend 
analytically across the boundary. 

7. The Subadditivity Problem for Analytic Capacity 

This section is about the study of the following question : 
Is it true that 

(6) j(E U F) < j(E) + 7 (F) 

for all compact sets F, F? 

Suita [16] proved that ([5]) holds if E, F are disjoint connected compact sets. 

One of the main obstacles in the study of inequality © is that it is difficult in 
practice to determine the analytic capacity of a given compact set. However, the 
numerical examples in the last section show that our method is very efficient when 
the compact in question is a finite union of disjoint disks. Fortunately, this partic- 
ular case is sufficient : 

Theorem 7.1. The following are equivalent : 

(1) j(E U F) < 7(F) + 7(F) for all compact sets F, F. 

(2) 7(F U F) < 7(F) + 7(F) for all disjoint compact sets F, F that are finite 
unions of disjoint closed disks, all with the same radius. 

Clearly (1) implies (2), but the fact that the converse holds is nontrivial. For the 
proof of Theorem 17. 1[ we need a discrete approach to analytic capacity introduced 
by Melnikov [9]. 
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7.1. Melnikov's discrete approach to Analytic Capacity. Let Zi, Z2, ■ • • , z n £ 

C and let ri,r2, . . . , r n be positive real numbers. Define Z :— (zi, Z2, ■ ■ ■ , z n ) and 
R := (ri, r%, . . . , r n ). Suppose in addition that \zj — Zk\ > Tj + Tk for j ^ k, so that 
the closed disks D(zj,rj) are pairwise disjoint. Set 

71 

K(Z,R) := \jB( Zj , rj ) 



and let 



/ii = nx{Z,R) := sup 



£< 



where the supremum is taken over all points ai, . . . , a„ € C such that 



< 1 



Clearly, we have [i\ < 7(K(Z, R)). 



Finally, for any compact set K C C and 6 > 0, we write Kg for the closed 
5— neighborhood of K. 

The following lemma is precisely what we need to prove Theorem 17.11 : 



Lemma 7.2. Let K C C compact, and let 5, e > 0. XTien i/iere exist zi, ■ ■ ■ , z n £ 
and < r < S such that \zj — Zk\ > 2r /or j 7^ fc, and 

m (Z,i?)>(l-e) 7 (X) 

where Z — (z\, . . . , z n ) and R — (r, . . . , r). In particular, 

j(K(Z,R)) > (l-e)j(K). 

Proof. See [HI Lemma 1]. 

We can now prove the reverse implication in Theorem 17.11 : 
Proof. Suppose that (2) holds, but that there exist compact sets E,F with 

j(EUF) > 1 {E)+ 1 {F). 
Let < e < j(E U F) - j(E) - 7(F). Take S > sufficiently small so that 



□ 



(7) 7(^25) <l(E)+e/3, 
and 

(8) 7(^25) < 7(F) + e/3. 

By Lemma 1731 there exist z\, z%, . . . , z n £ (EU F)g and < r < 6 such that 

1 (K(Z,R))> 1 (EUF)-e/3 

and the disks D(zj, r) are pairwise disjoint. For each j £ {1, 2, . . . , n}, fix w, € SUi* 1 
with \zj — Wj\ = dist(zj, 2? U F) < 6. Let A be the union of the disks U>(zj,r) with 
Wj £ E, and let B be the union of the disks D(zfc,r) with G F\E. Then 
A C _E 2(5 and B C F 2 s- Since 

7(iUB) > 7(£Uf)-e/3, 
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we have 

7(^4 U B) > 7 (EUF)- £ /3 

> 7(F) + 7(F) + e-e/3 
= 7(F) + e/3 + 7(F) + e/3 

> 7 (A)+ 7 (B), 

where we used equations ([7]) and ([5]). This contradicts (2). 



□ 



7.2. Discrete Analytic Capacity. For Z = (zx, . . . , z n ) G C™, where Zj ^ Zk for 
all j k, and r > 0, define 



j(Z,r) := 7 



U =1 



Assume in addition that the discs V>(Zj, r) are pairwise disjoint. 

By Theorem 17. 1[ the subadditivity of analytic capacity is equivalent to 

1 (Z,r)< 1 (Z',r)+ 1 (Z" 1 r) 

for all z\, . . . , z n € C and all m € {1, 2, . . . , n — 1}, where Z — (zi, . . . , z„), Z' = 
(zi,...,z m ) and Z" = (z m+ i, . . . , z n ). 
The above inequality can be written as 

R(Z,r,m) < 1, 

where 

R(Z,r,m) := — "4^-7 r- 

{ ' l{Z',r)+ 1 {Z" 1 r) 

The above shows the importance of studying the quantity R(Z, r, m). 

In this subsection, our objective is to obtain the following asymptotic expression 

for R(Z, r, m): 

Theorem 7.3. Fix z\, Z2, ■ ■ ■ , z n £ C and fix m G {1, 2, . . . ,n — 1}. Then 

R(Z,r,m) = 1 - Cr 2 +0(r 3 ) 

as r — > 0. where C is a strictly positive constant depending only on m, n and 
Zi, Zi-, ■ ■ ■ , Z n . 

For the proof of Theorem 17.31 we need to introduce a discrete version of analytic 
capacity, first considered by Melnikov [9]. 

For Z — {zi, . . . , z n ) and r > 0, define the discrete analytic capacity A by 

2 ' 



X(Z, r) := sup 



i=i 



where the supremum is taken over all (ax, . . . , a n ) S C n such that 



E 

k=l 



— 



< 1. 
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We also introduce the following constants: 



M{Z,r) 
N(Z,r) 



4-4- \z k -Zj\* 



1/2 



M(Z, i?) 



1/2 



k j^k 



Then we have 



Theorem 7.4. Let Z = (zi, . . . ,z n ) G C n and r > 0. Suppose that the discs 
V)(zj,2r) are pairwise disjoint. Then 

I J| f <A(Z 1 r)<(l + ^,r))7(Z 1 r). 

Corollary 7.5. Lei if C C compact, and let 5 > and e > 0. XTien i/iere exists 
Zi, . . . , z n G if 5 and < r < 5 such that \zj — Zk\ > 2r /or j ^ k and 

| 7 (A-)-7(2,r)| <e, 

where Z = (z±, . . . , z n ). Furthermore, Z and r can be chosen so that M(Z, r) < e 
and N{Z,r) < e. 

See Theorem 2]. 

Now, we need another expression for \(Z, r) which is easier to manipulate. We 
proceed as in [3] : 



It is not hard to show that 
» /■ 

E 



V 



E 



0/~ Zo 



2 \ 



/ 



((D^+B)a,a), 



where a = (01, . . . ,a n ), -D_r is the diagonal matrix with each entry of the main 
diagonal equal to r, and B = (bjk), where 

r 



h jk= E 



[Zj Zm)(zk Ztri) 



Here (■, ■) is the standard scalar product in C n . 

The matrix B can be written in the form B = CDrC*, where C — {cjk) is the 
Cauchy matrix associated with z%, . . . , z n , i.e. 

c jk = - (j ^ k), Cjj = 0. 

Zj Zfc 

Arguing as in []J[ Lemma 3], we obtain 

\{Z,r) = {{D R l +CD R C*)- 1 {\),1), 
where 1 = (1, 1, . . . , 1) G C". 



The following lemma contains estimates for the discrete analytic capacity: 
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Lemma 7.6. Let Z = (zi, . . . , z n ) 6 C™ and let r > 0. Then 

nr -r 3 (CC*{l),l) < \{Z,r) < nr - r 3 (CC* (1), 1) + r 5 (CC*CC* (1), 1). 

Proof. 9, Lemma 4]. 
We shall also need the following lemma : 

Lemma 7.7. Fix Zi,...,z n € C and, as before, let Z — (zi,...,z n ), Z' = 
(zi, ■ ■ .,z m ) and Z" = (z m+1 , . . .,z n ). 

Write a := (CC*{1), 1), a 1 (C'C'*{1), 1) and a" (C"C"*{1), 1) , where C is 
the Cauchy matrix associated with Z' and C" , the Cauchy matrix associated with 
Z". 

Then 

a > a' + a". 

Proof. First note that we have the following expression for a : 
(CC*(1),1) = (C*(1),C*(1)) 

= EEE 



E E u. _ r ,i2 + E E E 



An elementary calculation shows that the second sum is equal to 

\ ^ / 4S(zj,Zk,zi) \ 2 _ m—^ 1 
ji^<i \ \ z i ~ z k\\ z j ~zi\\z k -zi\) ,^ <t fl(zj-,2fc,z,) 2 ' 

where S(zj, z{) is the area of the triangle with vertices Zj, Zk, Z\ [S{zj, Zk,z{) = 
if Zj, Zk, zi are collinear) and R(zj, Zk,zi) is the radius of the circle through Zj, z k 
and z\ [R(zj, Zk, z{) = oo if Zj, Zk, Z\ are collinear). 

The conclusion follows. 

□ 

We can now proceed to the proof of Theorem [731 : 

Proof. Write M(Z, r) := Ar A and N(Z, r) := Br 3 , where A, B do not depend on r. 
Also, define 

:= (CC*CC*(1), 1), p := (C'C'*C'C'*(1), 1), P" := (C" C"* C" 'C* (1) , 1), 
where C, C , C" are as in Lemma 17771 
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By Theorem 17.41 and Lemma [721 we have 

(1 + 4Br 3 )(l + 2Ar 4 )\(Z, r) 



R(Z, r, m) < 



X(Z',r)+X(Z",r) 
< (1 + 4Br 3 ){l + 2Ar 4 ) 



, j nr — ar + (5r° 



mr — a'r 3 + (n — m)r — a"r 3 



„ , w ... n — ar 2 + Br 

{l + ABr 3 ){l + 2Ar 4 )- 



n — (a' + a")r 2 

Now, by Lemma [7771 a = a' + a" + S for some 5 > 0, and thus 



3\m . 4 ,n - (a' + a" + 5)r 2 + f3r 4 



R(Z,r,m) < (1 +450(1 + 2Ar 4 ) 



(1 + 45r 3 )(l + 2,4r 4 ) 1 - r 



r 4 



n — (a' + a")r 2 n — (a' + a")r 2 



= (l+4Br 3 )(l + 2Ar 4 ) (l--r 2 + 0(r 4 ) 

\ n 

= 1- -r 2 + (3(r 3 ), 
n 

as r —> 0. 

Proceeding similarly, we obtain the reverse inequality : 

4n-ih , ,o„3s-i A(Z,r) 



R(Z,r,m) > (l + 2Ar 4 )-\l+ABr 



X(Z',r)+X(Z",r) 

> (l + 2Ar 4 )" 1 (l+4 J Br 3 )~ 1 



mr — a'r 3 + f3'r 5 + (n — m)r — a!'r 3 + /3"r 5 

(i + 2Ar 4 r ( i + ^r \_ {a , +a :- r r + \ p , +nri 

. as t , t*. i n - (a' + a")r 2 — Sr 2 

(l + 2Ar 4 )~ 1 (l+4 J Br 3 ^ 1 



n-(a? + a")r 2 + (8' + /3")r 4 



n - (a' + a")r 2 + (8' + B")r 4 
(1 + 2Ar 4 )" 1 (l + 4Br 3 )" 1 ( 1 - -r 2 + 0(r 



n 
n 

as r — > 0. 

□ 



Remark 3. One immediately deduces from Theorem 17.31 some interesting prop- 
erties of the ratio R(Z,r,m). First of all, for fixed Z and m, the asymptotic 
expression implies that R{Z, r, ni) < 1 for all r sufficiently small. Second, we have 
that R(Z, r, m) — > 1 as r — > 0. Another interesting consequence of Theorem 17.31 
will be discussed in the next section. 
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8. A Conjecture related to the Subadditivity Problem 

8.1. Formulation of the conjecture. Recall that, by Theorem 17. 11 the subaddi- 
tivity of analytic capacity is equivalent to the subadditivity in the case of disjoint 
finite unions of disjoint disks, all with the same radius. For such compact sets, 
our method for the computation of analytic capacity is very efficient, see e.g. the 
numerical examples of Section [6l This allowed us to perform a lot of numerical 
experiments. 

More precisely, let Z = (zx,...,z n ) E C™ and m g {l,...,n— 1}. Let r > 
sufficiently small so that the disks B(zj,r) are pairwise disjoint. Recall that 
R{Z, r, to) is defined by 

where 

m 

E:=\jB( Zj ,r) 
j'=i 

and 

n 

F:= |J B( Zj ,r). 

j=m-\-l 

Using our numerical method, we can easily compute upper and lowers bounds for 
the ratio R{Z, r, m). 

All the numerical experiments that we did seem to suggest that analytic capacity 
is indeed subadditive, i.e. 

R(Z, r, m) < 1 

for all Z , r, m. More surprising though, all of our numerical experiments seem to 
indicate that the ratio R(Z, r, m) decreases as r increases. We formulate this as a 
conjecture : 

Conjecture. Fix Z = (zi, . . . , z n ) and m G {1, . . . , n — 1}. Then R(Z,r,m) is a 
decreasing function of r. 

In view of Theorem 17. 31 a proof of the above conjecture would imply that analytic 
capacity is subadditive. 

Moreover, Theorem 17.31 also implies that, for fixed Z and m, the above holds for 
all r sufficiently small. 

We now present several numerical experiments to illustrate the conjecture. 
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8.2. Numerical examples. 

Example 8.2.1. For the first example, the number of disks is n — 40, and m = 20. 
The compact set E is composed of the 20 disks with bold boundaries, and F is the 
union of the remaining disks. 



•oooooooo 
•oooooooo 
oooooooo 
oooooooo 

► oooooooo 



FIGURE 7. The set E U F for Example [8X11 



For 500 values of the radius r equally distributed between and 0.499, we computed 
lower and upper bounds for the ratio R(z,r,m). 

Below is the graph of the lower bound versus r. The graph for the upper bound is 
almost identical; the two graphs differ by at most 0.002481. 




0.1 0.2 0.3 0.4 0.5 

radius r 



FIGURE 8. Graph of the ratio j(E U F)/(i(E) + j(F)) for Example EJ 
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Example 8.2.2. Here the disks are centered at equally spaced points in the real 
axis. The number of disks is n = 10, and m = 5. 



FIGURE 9. The set E U F for Example 18X21 




0.1 0.2 0.3 0.4 0.5 

radius r 



FIGURE 10. Graph of the ratio 7 (£U F)/(-y(E) + 7(F)) for Ex- 
ample 
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Example 8.2.3. Here the disks are centered randomly. The number of disks is 
n = 18, and m = 12. 



7 $ 

o 



o 



o o 



o 



o 



o 
-0- 



o 



o 



o 



Figure 11. The set E U F for Example EZM 




0.1 0.2 0.3 0.4 0.5 

radius r 



FIGURE 12. Graph of the ratio 7 (£U F)/(-y(E) + 7(F)) for Ex- 
ample EH 
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Our last example shows that the situation can be different if, instead, we fix the 
radius of some of the disks and let the radius of the other disks vary. 

Example 8.2.4. Here the disks on the left are fixed, with radius 0.49. The radius 
of the two disks centered at 10 and 11 vary simultaneously, from to 0.499. 




Figure 13. The set E U F for Example [8^41 




FIGURE 14. Graph of the ratio ^{E U F)/(-y(E) + 7(F)) for Ex- 
ample IHHH 
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8.3. Proof of the conjecture in the case n = 2. We end this section by giving 
a proof of the conjecture in the simplest case. 

Theorem 8.1. Let E and F be disjoint closed disks of radius r. Then 

j{EUF) 
7 (£) + 7 (F) 

is a decreasing function of r . 

Proof. The main ideas of the proof that follows were suggested to us by Juan Arias 
de Reyna, and we gratefully acknowledge his contribution. 

Without loss of generality, we can suppose that E and F are centered at c and — c 
respectively, where c > and < r < c. In this case, we have, by formula ([3]) of 
Example 16.1.11 

l{EUF)= r -(±-^jV 2 {q)\ 
where q is the solution in (0, 1) of the equation 

Recall that the Jacobi theta-functions are defined by 



mq) ■■= e « (n+i/2)2 = v /4 n (! - 9 2n )(i + i 2n r, 

neZ n=l 

Mq) :=E9" 2 = na-'? 2 ")(i+^ 1 ) 2 J 

n£jj n—1 

oo 



for q G (0, 1) 



Now, since 7(£?) = 7(F) = r, we have that 
7(EUF) 1/1 



7( J B)+ 7 (F) 4 Vx/g / 

By equation ([9]). g increases from to 1 as r increases from to c. It thus suffices 
to show that f(q), defined above, is a decreasing function of q G (0,1). Proving 
this directly seems difficult, mainly because of the difference in the behavior of / 
near 1, and away from 1. For this reason, we separate the proof in two cases : 

Case 1 : q G (0, 0.8] 

The idea in this case is to express / as a Jacobi product and then compute the 
logarithmic derivative. Using the product expression for $2 (?) > we obtain 
00 00 
/(«) = (1 - «) A (1 - Q 2n f (1 + 9 2n ) 4 = (1 - g) n (1 - g 4 ") 2 (l + q 2n f 
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Since / is positive on (0, 1), proving that / is decreasing in this interval is equivalent 
to proving that its logarithmic derivative is negative. Computing the logarithmic 
derivative and multiplying by q, we obtain the function 

/'(g) Q ^ 8ng 4 " ^ Anq 2n 

u[q): ~ q m~ i- q h l -« 4n k l +i 2n ' 

We now estimate u. Fix k > 2. We have 



fc-i 



-W S -T^ - E T^r + E ^ - E W» + E 4n f 

n=l y n=l n=fc n=fc 

by splitting the sum between n < k — 1 and n> k, and using the inequalities 



2 71 



8n<7 4 " 

and 



> 8nq An 



< 4nq z 



l +q 2n 

Evaluating the two infinite series, we obtain the following upper bound for u(q): 



q 5^-4 8nq 4n ^4nq 2n 8(k - l)g 4 ( fe+1 ) - 8kq 4k Akq 2k - A{k - l) g 2 ( fe + 1 ) 



q \ ^ oil' 



+ - 



q rf=l 1 ~~ ^ n^l 1 + t 1 ~~ ^ ^ ~~ 

Using maple, we can substitute different values of fc and solve where the resulting 
expression is negative. With k = 10, we obtain that the above expression is nega- 
tive for q < 0.81121. In particular, / is decreasing in the interval (0, 0.8]. 

Case 2 : q e (0.8, 1) 

In this case, we shall make another change of variable, using the modularity of the 
theta-functions. The Jacobi modular identity for theta-functions implies that 

Making the change of variable q = e"^l x , we get 

/(g) = -^q- 1 ' 2 q^ 2 )Mq) 2 = \x (sinh |-) ^(e~™ ) 2 . 

Note that x increases from to oo as q increases from to 1. Furthermore, if 
0.8 < q < 1, then x > n/ (log 5/4) w 14.0788. It thus suffices to prove that f(x) is 
decreasing for x > 14. Write f(x) = g(x)h(x), where 

a ( x ) t:X sinh — 
yw 2 2x 

and 

h{x) := ^(e - ™) 2 . 

It is easy to prove that g is decreasing on (0, oo). Indeed, g'{x) < is equivalent to 

7T 7T 

tanh — < — , 
2x ~ 2x 

which is true since tanh(0) < 9 for 6 > 0. 

Now, we shall use the fact that h(x) w 1 for x > 14 to deduce that in this case, 
the behavior of / and g are nearly the same. We have to do some numerical error 
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analysis : 

First, we need to estimate how close h(x) is to 1 when x > 14. Note that 
t? 4 (g) = l-2q + 2q 4 -2q 9 + ... > 1 - 2q - 2q 2 - 2q 3 - . . . 

= 1 2(1 



Hence, for small q, 

M<if > i + ^rr^ 4-i- > i - 4-?-. 

{l-q) 2 l-q l-q 

Thus, with q = e~ nx , 

p-TTX 

1 - Mqf = 1 - t?4(e- TO ) 2 < 4- — < Be"™, 

1 — e 

since 1 - e" 7 ™ > 1/2 for ir > 14. Define £(x) := 1 - i^e" 7 ™) 2 , so that 

(10) < £(x) < 8e~ nx , 

for x > 14. 

We shall also need an estimate for the derivative of g : 

... I/-, 71 " n , \ 1 i 71 "/ i ^ ^ \ 

— g (x) = — sinh cosh — = - cosh — — tanh 1 . 

y w 2 V 2x 2x 2x/ 2 2x V 2x 2x/ 

Now, since x > 14, we have 

- cosh — < - cosh — w 0.50315 < 1. 
2 2x ~ 2 28 

Also, for > 0, 9 — tanh# < 9 3 /3. Indeed, both functions are at and if we 
compare the derivatives, we get 

tanh 2 (9 < 9 2 , 

which holds for every 9 > 0. 

We thus obtain the following estimate for the derivative of g : 

1 / 7T \ ^ ir 3 IT 3 

Now we estimate g(x). This is easy; g is decreasing, so for x > 14, we have 

(12) g(x) < g(U) = 1(14) sinh ^ « 0.78705 < 1. 
We also have the easy estimate 

(13) < !• 
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Finally, we estimate the derivative of $4 : 

00 

= 2(-l + 4< z 3 -9< Z 8 + 16< ? 15 -...) > 2(-l-^(n + l)g n ) 



n=l 

-2-4—^ + 2 q2 



> -2-4 



{l-qf {l-qf 

q 



Hence, 

(14) -tfAe-™) < 2 + 4^ < 2 + l&e-™ < 3, 

again for x > 14. 

We now have everything we need to estimate f'(x) : 
f'(x) = g'(x)h(x)+g(x)h'(x) 

= g'(x) - g'(x)£(x) - 2vg{x)Me-™W^-™)z-™- 
Hence, by equations ((TU)) to (fT4l : 

/'(as) < <?'(a;)+ 86^ + 2^(1) (1) (SK™ = ^(i) + 8e~™ + SkiT™ 

< g'(x) + 27e-™, 

and this holds for all x > 14. To complete the proof, all that remains is to prove 
that this last expression is negative for x > 14. We proceed as follows : 

g'(x) + 27e~ 7rx = - cosh — ftanh — - — ) + 27e~™ 
y w 2 2x V 2x 2x1 

< -ftanh )+27e- ,ra . 

~ 2 V 2x 2x) 

However, we have the inequality 

e 3 

(15) tanh6»<6» + 9° (9 > 0), 

3 

which follows from Taylor's theorem. Indeed, first note that 

d 5 (tanhy) 42 6 42 
— = — 88sech ytanh r/ + 16sech y + 16tanh ysech y. 

Now, since tanhy < 1 and sechy < 1, the remainder in Taylor's theorem for tanh# 
is less than (320 5 )/5! < # 5 , and inequality (fT5|) follows. 
With 9 = tt/(2x), we obtain 

. ^ 1 / 7T 1 / 7T \ 3 /7T\ 5 7T \ 

48 64 x 5 

Clearly the last expression is negative for x > 14, and thus we are done. 

□ 
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